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Abstract 



The recently reanalysed LSND data is investigated in the context of 
three flavor oscillations, with an emphasis on mass hierarchies and S13. 
The resultant mass hierarchies and oscillation angles are tested with 
Q_i| regard to "key" features of solar neutrinos and atmospheric neutrinos, 

O ■ e.g., "average" survival probability in the case of solar neutrinos, and 

zenith angle dependence and up-down asymmetry in the case of high 
energy atmospheric neutrinos. We find there are three distinct mass 
hierarchies, e.g., Am « AM, Am < AM and Am ~ AM. In the first 
^ . and second case, the calculated range of S13 is in agreement with the 

"LMA" solution of Akhmedov et al, the lower limit on S13 in these cases 
is also in agreement with the recent analysis of Garcia et al. based on the 
constraints of SNP, ANP and CHOOZ, therefore, strongly supporting 
the neutrino oscillations observed at LSND. Further, the solutions of 
S13 found in the third case correspond to the value of S13 found by 
Akhmedov et al. in the case of "SMA" and "LOW" solutions. A rough 
estimate of the possibility of the existence of CP violation in the leptonic 
sector is also carried out for different possible ranges of S13, indicating 
that the CP asymmetries may be measureable even in the case of LSND. 

Neutrino oscillations provide the preferred solution for the anomalous be- 
haviour observed in the case of solar neutrinos [0-0, atmospheric neutrinos 



||-[14j], as well as at LSND [ l5| . This essentially implies that the neutrinos 
are massive particles and the observed flavor eigenstates are linear combina- 
tions of mass eigenstates, in parallel to the quark mixing phenomenon. By 
far, the strongest evidence for neutrino oscillations is provided by the Su- 
perKamiokande(SK) atmospheric neutrino data ||. In view of the fact that 
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the signal for neutrino oscillations is also the signal for physics beyond the 
Standard Model (SM) |I§ |T7|, an intense amount of activity 



both at the experimental as well as at the phenomenological level, is going on 
to fix the mass hierarchies as well as the parameters of the neutrino oscillations. 

Several detailed and exhaustive analyses have been carried out in the 
case of the solar neutrino data as well as in the case of the atmospheric neutrino 
data. These analyses along with the results from several other experiments 
have provided valuable information about the masses as well as about the 
mixing parameters. The constraints on masses and mixings are presented in 
terms of the mixing angles (#12, #13, #23), defined analogous to the quark mixing 
angles p4| , and the mass square differences. The most important constraint 
is provided by the SK analysis of their atmospheric neutrino data, || ^5], , 
for example, 

Am 2 32 ~ (2 - 6) x W~ 3 eV 2 , sm6 23 ~ 0.545 - 0.839, (1) 

where subscripts '2' and '3' correspond to and v T . Although the constraint 
is for two flavor oscillations, however, in view of the CHOOZ constraint ]27[ 



sin^ia = \U e3 \ 2 < (0.06 - 0.018) for Am 2 31 = (1.5 - 5) x l0~ 6 eV% (2) 

it remains valid in the case of three flavor oscillations also. It is generally be- 
lieved that the bimaximal mixing matrix |28| provides a good approximation to 
the neutrino mixing matrix and is also consistent with the recent most exhaus- 
tive three flavor analysis of Solar Neutrino Problem (SNP) and Atmospheric 
Neutrino Problem (ANP) by Garcia et al. |26 . 



It is widely believed that the oscillation solutions to SNP, ANP and LSND 
data requires three different mass scales [p9fl , therefore a simultaneous solution 
seems to be possible only within four flavor(3 active + 1 sterile) scenario [25 



However, the latest analysis of the SK data || strongly indicates the absence of 
the role of sterile neutrinos in understanding SNP and ANP. This has compli- 
cated the issues related to neutrino mass scales. Besides this, there are several 
issues where more information is needed to deepen the understanding of the 
phenomenology of neutrino oscillations. One of the important issue is to find 
the lower limit of S13, the upper limit, of course, is provided by the CHOOZ 
constraint. Some attempts have been made in this direction |26], however 
a clear picture is yet to emerge. Similarly, despite considerable progress in 
understanding the neutrino mass differences, more information is very much 
desirable in knowing the various mass hierarchies which can fit the solar neu- 
trino data, the atmospheric neutrino data and the LSND data simultaneously. 
The other burning issue, in the context of neutrino oscillations, is the exis- 
tence or non existence of CP violation in the leptonic sector, for which serious 
effort is going on to decipher it in the Long BaseLine(LBL) experiments. Any 
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clue about CP violation, from phenomenological considerations, would have 
far reaching consequences for the neutrino oscillation phenomena as well as for 



the planned LBL experiments and Neutrino Factories |2(J, |21], |22], pi ]. 

The reaffirmation of the LSND results presented at Neutrino 2000 as well 
as the lowering of their probability, would have important consequences as 
this is an appearance experiment, unlike solar neutrino experiments and at- 
mospheric neutrino experiments. Another extremely important observation 
which has been emphasized at Neutrino 2000 is the absence of energy depen- 
dence of solar neutrino recoil energy spectrum by SK || as well as by Sudbery 
Neutrino Observatory (SNO) [[jj. Similarly, the preliminary observations of 
K2K experiment |32| as well as of SNO endorse the conclusions of SK regard- 
ing the hypothesis of neutrino oscillations. In view of the recent observations 
of SK about sterile neutrinos as well as the issue of three different mass scales 
being required to explain SNP, ANP and LSND results, the question regarding 
the number of neutrino species has become extremely important. In this con- 
text, Minakata has strongly advocated that one should try to re-explore, 
within three active flavors, the crucial features of the solar neutrino data, the 
atmospheric neutrino data and the LSND data which can be reconciled si- 
multaneously. Some preliminary attempts |33], [34], ^] have been made in this 
connection which need to be updated in the context of the observations at 
Neutrino 2000. Further, these attempts have not gone into the details of the 
mixing angles, in particular the lower limit of S13 and its implications on CP 
violation. Furthermore, these analyses do not adequately emphasize the is- 
sues related to flux uncertainties as well as up-down asymmetry in the case of 
atmospheric neutrinos. 

The purpose of the present communication, on the one hand, is to explore, 
within three active flavors, the possibility of reconciling the latest LSND data 
with some of the "key" features of solar and atmospheric neutrino data. On 
the other hand, it is to explore the constraints on the element S13, the resultant 
mass hierarchies as well as the possible clues about the order of CP violation 
in the leptonic sector. 

Before we proceed further, we underline the "key" features of the solar 
neutrinos and the atmospheric neutrinos which are relevant to the present 
investigation. In the case of solar neutrinos, the most important feature of 
the data is the ratio of flux of solar neutrinos observed by various ongoing 
experiments to the flux predicted by Standard Solar Model(SSM) p5| . In 
view of the observation at SK and SNO about the energy independence of the 
solar neutrino flux, it is expected that the above mentioned ratio should be 
same in all the solar neutrino experiments. In this situation, as advocated by 
Bahcall et al. [36|, and Conforto et al. ||34|| , one can consider the "average" 
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solar survival probability as, 



P. 



sol 



0.50 ±0.06. 



(3) 



Similarly, in the case of atmospheric neutrinos the "key" aspects of the data 
which indicate the possibility of neutrino oscillations are zenith angle depen- 
dence of the ratio (observed versus Monte Carlo (MC)) of the neutrino fluxes 
and the up-down asymmetry. In view of the fact that the ratios are very much 
dependent on the MC simulations of neutrino fluxes, it is perhaps desirable 
to concentrate on those features of the data wherein the flux uncertainties are 
minimum. Details in this regard will be discussed later in the text. 

To begin with, we consider the neutrino mixing matrix, 
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U n U 12 U 13 
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U 31 U 32 C/33 




(4) 



where u e , u^, v T are the flavor eigenstates and Ui, i/ 2 , ^3 are the mass eigen- 



states. In the PDG representation | 24| , the mixing matrix can be expressed 
as . 



^ C12C13 
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(5) 



In the sequel, we detail the relevant appearance and disappearance proba- 
bilities in terms of masses and mixing angles pertaining to the LSND experi- 
ment, the SNP and the ANP. In the case of LSND experiment, we express the 
probability in terms of the mixing matrix. The expression for Plsnd, ignoring 



CP violation effects, can be expressed as p7 | 

-4t72it/iit/22£/i2sin 2 
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(6) 



L denotes 



where Am 2 = Am^ = m\ — ra\, AM 2 
the neutrino flight path i.e. the distance between the neutrino source and the 
dectector and E is the average energy of the neutrinos. If Am 2 , AM 2 are 
expressed in eV 2 and L is expressed in metres then E is in MeV. The LSND 
set up is characterised by L = 30m and 36MeV < E < 60MeV. For the 
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purpose of calculations, we take the average energy E = A2MeV and from 
their latest analysis |T5 |, we have 

P LSND = (2.5 ± 0.6 ± 0.4) x 10" 3 with Am 2 > 0.2eV 2 . (7) 

It is to be noted that when the above constraint is combined with the limit 
provided by the CDHSW data f38|, the range of Am 2 LSND becomes 

0.2eV 2 < Am 2 LSND < OAeV 2 . (8) 

In view of the fact that Su and S23 are fairly well known, one can make use 
of the equation with the identification of AM 2 with Am 2 LSND , to find the 
element Si 3 for different mass hierarchies, however these values have to satisfy 
the "key" features of the solar neutrinos as well as the atmospheric neutrinos. 

To this end, we first consider the case of solar neutrinos. The survival 
probability of solar electron neutrinos is given by, 

/ Am^ J , \ 

P sol = 1 - AUnUnUuUusm 2 ' " > 



V 4£ 

,2 



., ,(Am 2 + AM 2 )L\ 

-AU 12 U 12 U 13 U 13 sm 2 (^j^J ■ (9) 

The extremely large L/E factor, as well as the range of Am 2 considered here, 
enables one to simplify the above expression by carrying out averaging, which 
in terms of the mixing angles becomes 

p sol = 1-24(4 + 444). (10) 

In the case of atmospheric neutrinos, the most important feature is the 
zenith angle dependence of the ratio of the observed neutrinos to the MC 
expectations. For example, the ratios in the case of electron neutrinos (R e ) 
and muon neutrinos (R^) are defined as follows, 

R e = P ee + rP fie , (11) 

R^i — P^l + -Pfxe.1 (12) 

where P ee , P^ are the survival probabilities, P^ e is the transition probability, 
and 

N° 

r = TT ( 13 ) 

N? v ' 
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is the expected ratio of fluxes without oscillations (where N",^ is the initial 
flux of muon (electron) neutrinos). The ratio r varies from 1.6 ( for sub GeV 
neutrinos ) to 3.0 ( for multi GeV neutrinos ). Further, this ratio depends on 
the zenith angle also, varying from 1.6 for the horizontal neutrinos to 3.0 for 
the vertical ones. It is obvious that the uncertainties in r are bound to affect 
the calculation of the zenith angle dependence of R e and R^. Since the purpose 
of the present communication is not to go into the extensive details of these 
aspects of data, therefore we would like to concentrate on that part of the data 
which is free from such uncertainties. Interestingly, it is to be noted that these 
uncertainties get highly suppressed for the case of high energy upward going 
neutrinos. It is also to be noted that the SK detector depicts better correlation 
for upgoing high energy neutrinos with the charged leptons produced in the 
detector. Further, the asymmetry A, defined as 

where U is the number of upward going events with zenith angles in the range 
— 1 < cos# < —0.2 and D is the number of downward going events with 
0.2 < cos# < 1, is better defined for high energy neutrinos . 



In view of the above mentioned reasons, it is desirable to consider the zenith 
angle dependence of R e and R^ for high energy upgoing neutrinos. The zenith 
angle dependence is generally expressed in terms of L/E. It is to be noted that 
for the vertically upgoing neutrinos, the length (L) traversed is ~ 13000 km. 
For the SK geometry, the horizontal neutrinos correspond to those travelling ~ 
500 km, therefore, for the purpose of studying zenith angle dependence which 
does not involve uncertainty due to flux ratio r, it is desirable to avoid lengths 
~ 1000 km. 

Keeping in mind the mentioned constraints, the three flavor probability 
expressions P ee , ? w and P^ e required to evaluate R e and R^ can be expressed 
as follows. The electron neutrino survival probability is given by, 

P ee = l-4U u U u U 12 U 12 sm 21 ' A '"~M 



AE 



-AUnUnUuUuam 2 



(Am 2 + AM 2 )!/ 



AE 



-AU 12 U 12 U 13 U 13 sin 2 (^^j ■ (15) 
Similarly the muon neutrino survival probability is given by, 
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Pan = 1 - AU 21 U 21 U 22 U 22 sin 2 



-AU 2l U 21 U 23 U 23 sm 2 



Am 2 L\ 

{Am 2 + AM 2 )V 
4E 



-AU 22 U 22 U 23 U 23 sin 2 f^^J • (16) 



Likewise, the transition probability is expressed as, 

( Am 2 V 



P„ e = -4U 21 U n U 22 U 12 sm 2 



-AU 21 U n U 23 U 13 sm 



\ 4£ J 
2 /(Am 2 + AM 2 )V 



AE 



-4U 22 U l2 U 23 U 13 sm 2 (^j^J • (17) 

The factor L/E is quite large for the upgoing neutrinos considered here, there- 
fore, for the purpose of calculations the above expressions get reduced to the 
following, 

P ee = 1 - 2U 2 n U 2 2 - 2U 2 u U 2 l3 - 2U 2 l2 U 2 l3 , (18) 
P m = 1 - 2U 2 1 U 2 2 - 2U 2 21 U 2 2Z - 2Ul 2 Ul 3 , (19) 

P^ = 2U 2 3 U 2 3 - 2U 22 U 12 U 21 U n . (20) 

In Tables [TJ and ^, we have presented the results of our calculations, for dif- 
ferent mass hierarchies, incorporating the constraints due to the LSND results 
as well as due to the "key" features of SNP and ANP. To begin with, using 
the LSND probability given by the equation (|>|), we have calculated s\ 3 . The 
values of AM 2 are constrained by the equation (§), the values of s 23 by the 
equation (j^) whereas the values of su used for calculations are fairly consistent 
with bimaximal mixing as well as with the analysis of Garcia et al. |26 . 



In Table [I], we have presented a representative sample of the calculations 
carried out in the case of "natural" mass hierarchy Am 2 AM 2 , with Am 2 = 
10~ 4 eV 2 . The results corresponding to the hierarchy Am 2 < AM 2 are very 
much similar to those of Am 2 AM 2 , hence are not presented in the tables. It 
is interesting to mention that while carrying out the analysis of the range of Si 3 , 
using the LSND constraint, we find that s 12 plays no role. However, because 
of the fact that P so i depends on su values, it is not possible to consider values 
of S12 < 0.6. The fact that Plsnd is independent of s\ 2 values can be easily 
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understood if one carefully examines the equation (|]). Since Am 2 <C AM 2 , 
the first term of the equation (ffl) can be neglected and the contribution of Am 2 
from the second term can be ignored. The second and the third term can be 
combined, yielding in terms of mixing angles, 



D LSND = 4S^C 2 3 S 2 3 S 




(21) 



which has no S12 dependence. Therefore, while carrying out the calculations, 
we have considered the bimaximal mixing value of S12. From Table [1], it is quite 
evident that when we scan the possible range of Plsnd, AM 2 and s 23 , we get 
the range of S13 as (0.07 — 0.29). From the table it should be noted that P so i 
is very much in agreement with equation (||) for values of S13 only upto 0.18. 
Similarly, for the case of atmospheric neutrinos, R e and also agree with 
the data |n| for the same upper limit of S13. Therefore, when we consider the 



constraints imposed by SNP, ANP and the LSND experiment simultaneously, 
the range of Si 3 gets restricted to (0.07 — 0.18). It needs to be added that for 
the region of L/E considered here, R e and R^ assume constant values because 
of averaging of various oscillatory terms in the probability expressions. 

In Table 0, we have presented the results of the calculations corresponding 
to the "degenerate" case, i.e. where Am 2 and AM 2 are of the same order of 
magnitude. As a first step, we have tried to find the largest possible value of 
Am 2 which can fit the various constraints imposed by LSND, SNP and ANP. 
In this context, we find that a value of Am 2 > 0.085 eV 2 is not able to fit 
Plsnd- It is also to be noted that the range of S23 is very much limited to 
values ~ 0.82, which is almost near the largest value admitted by SK. Using 
these values of Am 2 and s 2 3, we have calculated possible values of Si 3 as well 
as of S12 which fit the data. As is evident from the table, the values of S12 are 
also quite limited to the range 0.66 — 0.70, whereas, S13 takes on much smaller 
values, 0.001 — 0.03, as compared to the previous case where Am 2 -C AM 2 . 

A closer scrutiny of our results reveals several points. It needs to be em- 
phasized that the lower value of S13 found in the case Am 2 <C AM 2 , is in 
agreement with the value found recently by Garcia et al. ||26|| , by carrying out 
a simultaneous analysis of solar, atmospheric and CHOOZ data. Interestingly, 
the range of S13 found by our calculations is in complete agreement with the 



range found by Akhmedov et al. |3(|, in the case of "LMA" solution of the 
SNP. Further, the upper range of Si 3 , surprisingly, is in agreement with that 
of CHOOZ constraint, however, the present upper limit, unlike the CHOOZ 
case, is scale independent. Similarly, in the second case where Am 2 ~ AM 2 , 
the calculated values of S13 correspond to the value of S13 found by Akhmedov 
et al. pCfl , in the case of "SMA" and "LOW" solutions. This agreement of 
s 13 range found primarily from LSND data with other diverse analyses lends 
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strong support to the neutrino oscillations observed at LSND. 

In the above mentioned tables, we have presented the results for two dif- 
ferent mass hierarchies, e.g., Am 2 <C AM 2 and Am 2 ~ AM 2 . In the latter 
case, we have already mentioned a specific value of Am 2 which in fact is the 
highest value fitting the LSND data. In the first case, we have given a typical 
value, however our results are very much valid if Am 2 corresponds to either 
the "LMA" or the "LOW" solutions to the SNP. The calculations also remain 
valid even if Am 2 ~ 10~ 2 eV 2 . 

In the tables, we have not included the results of the asymmetries calculated 
from equation fll4|) . In case we impose the range of A within la, e.g.,- 0.35 
to - 0.25 ||, we find that the results corresponding to the degenerate case are 
within this range, however in the case Am 2 <C AM 2 only values of S13 < 0.1 
satisfy the above mentioned range of the asymmetry. 

It is interesting to mention that the present calculations for both the mass 
hierarchies are fully compatible with the disappearance observations at 
K2K. Although the K2K results show up only at la level, still their compati- 
bility with the present calculations supports the oscillation hypothesis. 

In the present calculations we have not considered the downgoing atmo- 
spheric neutrinos, primarily because of the flux uncertainties associated with 
these neutrinos as well as because of the fact that high energy neutrinos do 
not oscillate for these lengths. However, it may be mentioned that for the low 
energy neutrinos Am 2 L/E can be ~ 1, in which case, the calculations become 
highly scale sensitive. 

It is interesting to explore which additional features as well as mass hier- 
archies, apart from the "key" features of SNP and ANP, would be describable 
if one extends the present 3 active flavors to 3 active + 1 sterile scenario. The 
details of these will be published elsewhere. 

After having found mixing angles which are compatible with the three neu- 
trino experiments, it is desirable to examine the kind of texture specific mass 
matrices which would reproduce these mixing angles, as well as the possibility 
of CP violation in the leptonic sector. Regarding the mass matrices, some pre- 



liminary investigations have been carried out earlier [ 3S| , details in the present 
context would be published elsewhere. Regarding CP violation, assuming that 
S is nonzero in the leptonic mixing matrix, one can estimate the order of CP vi- 



olation, by calculating Jarlskog's rephasing invariant parameter J ||40 |, through 
the relationship 

I J I = 2 x area of the unitarity triangle, (22) 



for details we refer the reader to reference [41]. In analogy to the quark sector 



J can be evaluated by using any of the six unitarity triangles in the leptonic 
sector, through equation (p2). A nonzero area of any of the unitarity triangle 
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would imply CP violation, which we ensure by putting the constraint that the 
sum of the two sides of the triangle is always greater than the third side. Using 
this constraint J can be evaluated for any of the six unitarity triangles, however 
we will present the results corresponding to the unitarity triangle formed by 
considering the orthogonality of the I and II row of the mixing matrix. In 
figure 0, we have plotted the probability distribution of J by varying s 12 , $23 
and S13 within the ranges mentioned in the I row of the Table ^ and 5 in 
the range 0° to 180°. Interestingly the histogram shows a sharp peak around 
J = 2.126 x 10~ 4 , a corresponding guassian fit results in the range 

J = (2.126 ± 0.272) x 10~ 4 . (23) 

The corresponding most probable range for 5 is 

S = 1° to 10° (in I quadrant), 

170° to 180° (in II quadrant). (24) 

We have repeated the above procedure for the mixing angles listed in row II 
and III of Table |[ The corresponding J and 5 are also listed in the table. 
In the case S13 = 0.001 — 0.10, we have mentioned for 5 the value 90° ± 14°, 
however, 5 can also exist with significant probability in the range 1° to 5° in 
the first quadrant and correspondingly 175° to 180° in the second quadrant. 
It is interesting to consider the CP asymmetry in the case AM 2 ~ Am 2 . 



The CP asymmetry |2l| is defined as follows 



P(u^^u e ) - P(u^V e ) 4sin 2 # 12 sin<5 . (2Am\ 2 L\ 

For the degenerate case, as well as using the LSND parameters L/E, we get a 
rough estimate of Aqp as follows, 



i-CP 



sin 5. (26) 



As shown in the Table ||], for this case S can take values 1° to 10°, therefore the 
asymmetry can be reasonable even for short baseline experiment like LSND, 
however in the case of LBL experiments the asymmetry can take large values. 
Similarly, in the case of natural mass hierarchy, the largest value of Aqp for 
LBL experiments, for example K2K experiments, can be given as 



0.4. (27) 



In conclusion, we would like to emphasize that we have carried out a si- 
multaneous analysis, with three flavor oscillations, of the recently updated 
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LSND data along with the "key" features of solar neutrinos and atmospheric 
neutrinos. Interestingly, we find that there are three mass hierarchies which 
are able to fit the LSND data, the "averaged" solar neutrino probability and 
the zenith angle dependence for high energy neutrinos. The mass hierarchies 
found are Am 2 <C AM 2 , Am 2 < AM 2 and Am 2 ~ AM 2 , the latter being 
valid only for values of S23 — 0.82. Out of these three mass hierarchies, the 
first and the second possibilities give similar results. In a particular case with 
Am 2 = 10~ 4 eV 2 , we find the range of S13 to be 0.07 — 0.18. This range is in 



fair agreement with the "LMA" solution of Akhmedov et al. |30[ . Further, the 
lower limit of this range is in agreement with the most recent and exhaustive 
analysis by Garcia et al. based on a simultaneous analysis of solar neutrino 
data, atmospheric neutrino data and CHOOZ results. Furthermore, the range 
of S13, for example 0.001 — 0.032, found in the degenerate case corresponds to 
the values of S13 found by Akhmedov et al. in the case of "SMA" and 
"LOW" solutions. These observations, in a way, show the compatibility of the 
oscillations observed at LSND with the solar and the atmospheric neutrino 
data. If in the case of atmospheric data, the up-down asymmetry is also in- 
cluded in our analysis, then the solution space gets restricted if the results of 
asymmetry are included upto la level. However, if the results are included 
upto 2<T level, then most of the solutions remain valid. We have also made 
a rough estimate of CP violating asymmetries assuming the existence of CP 
violation in the leptonic sector. Interestingly, CP violation may be observable 
even at LSND. 
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Table 1: Calculated values of Si 3 , P a tm, P S oU Re-, R/j. in terms of 
S12, s 23 , Am 2 (eV 2 ), AM 2 (eV 2 ), with s 12 = 0.70. 
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Table 2: Calculated values of Si 3 , P a tm, P S oi, Re, R^ in terms of 
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Table 3: J and 5 corresponding to different ranges of S13. 
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Figure 1: Histogram of |J| plotted by considering the orthogonality relation- 
ship of I and II row of the mixing matrix. 
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